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In this work we report the study of conformation dependent electronic transport properties of 
DNA double-helix within tight-binding framework including its helical symmetry. We have studied 
the changes in localization properties of DNA as we alter the number of stacked bases within a pitch 
of the double-helix keeping the total number of nucleotides in the DNA chain fixed. We take three 
DNA sequences, two of them are periodic and one is random and observe that localization length 
increases as we increase the radius of DNA double-helix i.e., number of nucleotides within a pitch. 

We have also investigated the effect of backbone energetic on the I-V response of the system and 
we find that in presence of helical symmetry, depending on the interplay of conformal variation and 
disorder strength DNA can be found in either metallic or semiconducting and even in an insulating 
phase, which in turn successfully explain all the experimental hndings by a single model. 

PACS numbers: 72.15.Rn, 73.23.-b, 73.63.-b, 87.14.gk 


I. INTRODUCTION 

The advancements of nanoscience and technology with 
everyday encouraging a growing number of scientists 
across the various disciplines to devise ingenious ways 
for decreasing the size and increasing the performance of 
the nano-electronic circuits. One of the promising route 
is to use molecules and molecular structures as a compo¬ 
nent of those circuits. From these efforts a new branch 
has emerged called molecular electronics. Among dif¬ 
ferent branches of molecular electronics, DNA and alike 
biomolecules have drawn maximum attention in the last 
decade from both the theoreticians as well as experimen¬ 
talists and still growing in numbers. The main reason 
behind this attraction is the potential of DNA to become 
an inevitable agent for the future nanoelectronic devices 
and computers, as it might can serve in different ways in 
a nano-electronic circuits such as a wire, transistor or a 
switch depending on its electronic properties [l|, 0. Not 
only this, a precise knowledge of charge transfer mecha¬ 
nism through DNA could help in understanding the pro¬ 
cess like oxidative damage sensing, protein binding, gene 
regulation and cell division. On the other hand electrical 
properties, specially conductivity of DNA can be used for 
marker-free gene test Q which is one of the most highly 
desired biophysical methods 0. Inspite of the vast efforts 
from physicists as well as biologists around the world, 
charge transport results through DNA are still quite con¬ 
troversial 1MI|. Experimentally it is found that DNA 
can behave either as a g ood conductor semiconduc¬ 
tor fa [Till , insulator [l^, [l3 and even as a superconduc¬ 
tor [l3|| at low temperature. Several exp eriments both 
on synthetic periodic DNA chains Ig, [ll| as well as un¬ 
ordered sequence of basepairs [l3, [3 show the presence 
of a conduction gap in I-V curves at room temperature. 
Whereas linear response observed in Ref. [l^ and both 
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the staircase and linear behaviour in I-V curves shown in 
poly(dG)-poly(dC) chains Due to this experimental 
ambiguity and lack of understanding of charge transfer 
mechanism in DNA, leads to different phenomenologi¬ 
cal models in which charge transfer is mediated via po- 
larons [l^, solitons M or electrons or holes @,123,1211. 

This diversity of experimental findings on transport 
properties of DNA is due to several reasons such as, 
DNA varies widely in terms of its composition, length 
and structure, presence of counterions and impurities 
which can attach to the phosphate group of the back¬ 
bones, environmental effects, thermal vibration and con¬ 
tact resistance variation. In this communication, we try 
to address the effects of structure of DNA i.e., conformal 
behaviour on its transport properties. Experiments done 
more than half a century ago by Wilkins et. ai, (2^ first 
suggested that overstretched DNA (quite longer than its 
natural length) undergoes transition to a structure that 
can accommodate elongation up to twice the length of 
a relaxed DNA. Crucial developments in understanding 
mechanical properties of DNA was achieved via stretch¬ 
ing experiments [23 - 1^ . Depending on the stretching 
force applied, DNA first uncoils, then exhibit stiff elas¬ 
tic response and at last undergoes an abrupt structural 
transformation. Now as all the DNA are twisted (nat¬ 
ural double-helix structure) and the amount of twist- 
stretch i.e., radius of the helix varies from one situation 
to another, this study has to be made in details. People 
have already tried to study the effects of conformation 
introducing twist angle or chirality [l^, [l^ into ab-initio 
calculations. Studies also have been done on electronic 
properties of stretched DNA [2^ but the effects of heli¬ 
cal structure and conformality on its transport properties 
are yet not well explored. While study within much sim¬ 
pler tight-binding framework is hardly available in cur¬ 
rent literature. In our work we try to find out these ef¬ 
fects within the tight-binding model. To do this we follow 
Ref. , where a mechanical model of DNA is proposed. 
DNA being modelled as an elastic rod, wrapped helically 
by a stiff wire. The radius of elastic rod can change upon 
stretching with a Poisson’s ratio 77=0.5. The outer wire is 
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affixed to the rod helically with a given pitch. As stretch¬ 
ing force being applied, the elastic rod elongates in the 
length and its radius decreases. As a result the stiff wire 
overwinds and the number of turn increases. We take this 
mechanical model and interpret in the language of tight- 
binding formulation. We use twisted ladder model [30j, 
to imitate this mechanical model which includes both 
the helical symmetry and conformation. We have been 
able to show three different phases of DNA i.e., metallic, 
semiconducting and insulating depending on helical sym¬ 
metry, conformation (twist-stretching) and disorder. We 
have also found some structural configurations at which 
system hardly disturbed by external changes. 

This paper is organized in the following way: In Sec. II 
we discuss about our theoretical formulation and describe 
the model Hamiltonian. We explain our numerical results 
in Sec. Ill and summarized in Sec. IV. 


II. MODEL AND THEORETICAL 
FORMULATION 

DNA, carrier of genetic code of all forms of life, a tt- 
stacked array of four different nitrogenous bases adenine 
(A), guanine (G), cytosine (C) and thymine (T) attached 
among themselves via hydrogen bond following comple¬ 
mentary base pairing and coupled with sugar-phosphate 
backbones forming the double-helix structure. In most 
of the theoretical models, electronic conduction [sTI - f^ 
is assumed through the long-axis of the DNA molecule. 
To model DNA, in our present study, we take the tight- 
binding (TB) dangling backbone ladder model 
and add extra hopping channels due to the proximity of 
bases in the upper strand with the corresponding bases 
of the lower strand in the next pitch to incorporate its 
helical symmetry. The Hamiltonian for the said model 
can be expressed as (for schematic representation of this 
model we refer to [^ 1 

HdNA — ^ladder T ^helicity T ^backbone i (1) 

where. 
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where c\j and Cy are the electron creation and annihi¬ 
lation operators at the ith nucleotide at the jth stand, 
tij = nearest neighbour hopping amplitude between nu¬ 
cleotides along the jth branch of the ladder, = on¬ 
site energy of the nucleotides, = on-site energy of 
the backbone site adjacent to ith nucleotide of the jth 
strand with q{j) =t,4. representing the upper and lower 
strands respectively, = hopping amplitude between 
a nucleotide and the corresponding backbone site, v = 
interstrand hopping integral between nucleotides in two 
strands of ladder within a given pitch, v' — interstrand 
hopping integral between neighboring atomic sites in the 
adjacent pitches which actually accounts for the helical 
structure of DNA. Here n denotes the number of sites in 
each strand within a given pitch. For simplicity, we set 
= Cfe, tij = U and = 4- 
To explore the transport properties of DNA, we use 
semi-infinite ID chains as source (S) and drain (D) elec¬ 
trodes connected to alternative strands of the DNA in 
cross-wise fashion to the left and right ends respectively 
and the Hamiltonian of the entire system is given by 
H = Hdna + Hs + Hd + Htun- The explicit form of 
iJs, Hd and Htun are 


Hs 

= E (eclci + teJ+iQ + H.C.) , 

2= — 00 

(5) 

Hd 

00 

= E (ec4i+te|-FiG + H.c.) , 

(6) 


2=A^+1 


Htun 

= r fcjci -1- cj^CAT+i -1- H.C.) , 

(7) 


where r is the tunneling matrix element between DNA 
and the electrodes. 

To obtain transmission probability T{E) of elec¬ 
trons m, through DNA double-helix for this two- 
probe set up, we use the Green’s function formalism. 
The single particle retarded Green’s function operator 
representing the complete system le., ds-DNA and two 
semi-infinite electrodes, at an energy E can be written 
as = {E — H + where 77 —>■ O'*" and H is the 

Hamiltonian of the entire system. Using Fisher-Lee [ 3 ^ 
0 relation the two terminal transmission probability is 
defined as T{E) = Tr[rLG’TflG“], where E being the 
incident electron energy and the trace is over the re¬ 
duced Hilbert space spanned by the DNA molecule. The 
effective Green’s functions can be expressed in the re¬ 
duced Hilbert space in terms of the self-energies of the 
source and drain electrodes G’’ = [G“]'l' = [E — Hdna — 

where £^^“2) = ^tun'^^i'^tun and 
^SiD) = retarded 

(advanced) Green’s function for the source (drain) elec¬ 
trodes. Here and are the retarded and ad¬ 

vanced self-energies of the source (drain) electrodes due 
to its coupling with the DNA molecule. It can easily be 
shown that the coupling matrices rs(£)) corresponding to 
the couplings of the DNA chain to the source (drain) elec- 
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FIG. 1: (Color online). Lyapunov exponent ( 7 ) vs v' for three DNA sequences at several disorder strengths (w), for four 
different values of n = 3, 5, 7, 10. 7 decreses with increasing values of n for all the sequences irrespective of disorder strength, 
though the features of localization curves are claerly distinguishable for different sequences. There is no distinct changes for 
the critical values of v' (say, u(,) which corresponds to the minima of 7 with n. 
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trodes ^s(d) = —2 Whereas the self-energies 

are the sum of S^(.^j=As(£,)-fiAs(D), ^siD) being the 
real part of corresponds to the shift of energy lev¬ 
els of DNA, and the imaginary part is liable for 

the broadening of these levels. 

Considering linear transport regime, at absolute zero 
temperature, the two terminal Landauer conductance is 
given hy g = ^T{Ep), and the current passing through 
the DNA chain for an applied bias voltage V can be writ¬ 
ten as 


I{V) = - / TiE)dE , (8) 

“ JEF-eV/2 

where Ep being the Fermi energy. Here we have as¬ 
sumed that entire voltage drop occurs only at the bound¬ 
aries of the conductor. 


III. RESULTS AND DISCUSSIONS 





FIG. 2: (Color online). Lyapunov exponent ( 7 ) vs disorder 
(w) for three DNA sequences with i;'=0.3 eV, for four different 
values of n = 3, 5, 7, 10. Unifrom behaviour of localization 
has been observed for all the sequences for whole range of 
disorder. 

We first study the localization properties of the system 
by altering the number of bases in a given pitch of the he¬ 
lical structure. In order to do that we define localization 
length (1) from Lyapunov exponent (7) (^ . 

7 = 1/1 = — lim i < ln(r(£l)) > , (9) 

L—¥00 Li 

where L = length of the entire DNA chain in terms of 
basepairs, and <> denotes average over different dis¬ 
order configurations. Though other distribution func¬ 
tions e.g., Gaussian and binary have been used to simu¬ 
late experimental effects in previous studies [s^, but we 
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FIG. 3: (Color online). Lyapunov exponent ( 7 ) vs number 
of nucleotides within a pitch (n) for three DNA sequences 
with u'^0.3 eV, at different disorder strengths (w). Variation 
is quite uniform except for poly(dA)-poly(dT) sequence, a 
sharp peak is present there around n = 5 for heigher values of 
disorder, showing this may be the most localized configuration 
for that sequence. 


think it is appropriate to employ the most disordered case 
to simulate the actual experimental complications where 
the on-site energies of backbones Cf, to be randomly dis¬ 
tributed within the range [eb-w/2, eb-|-w/2], where Ih is 
the average backbone site energy and w represents the 
backbone disorder strength. For the purpose of numer¬ 
ical investigation the on-site energies of the nucleotides 
are chosen as the ionization potentials of the respective 
bases, i.e., ea = —0.56eU, ca = —O.OTeU, ec = 0.56eU, 
ct = 0.83eU. The intrastrand hopping integrals be¬ 
tween identical nucleotides are taken as t = 0.35eU 
while those between different nucleotides are taken as 
t = 0.17eU. We take interstrand hopping parameter to 
be 77 = 0.3eU. We emphasize that in case of the extended 
ladder model [12,113) diagonal hopping between different 
nucleotides are also taken into account. But as in our case 
no diagonal hopping being considered, we compensate 
this by taking a quite larger value of interstrand hopping 
parameter v. Now as all the nucleotides are connected 
with sugar-phosphate backbones by identical C-N bonds, 
we take the hopping parameter between a base and corre¬ 
sponding backbone site same for all 4 = 0.7eU 34] • The 
parameters used here are the same as those used in 14411 
which are consistent with ab initio calculations |45l - l47j| . 
For interstrand hopping v' between nucleotides of adja¬ 
cent pitches we follow Ref [13. Nevertheless, we want 
to mention that choice of the tight-binding parameters 
is not unique and several parameter sets have been pro¬ 
posed in the existing literature [i^ . 

In Fig. [1] we have plotted the variation of inverse lo¬ 
calization length (7) for three sequences with v' (which 
accounts for the helicity of DNA) at different values of n 
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FIG. 4: (Color online). Variation of Lyapunov exponent 7 with energy (E) for three DNA sequences with w'=0.3 eV, for 
different values of n. Effect of conformation (n) is stronger at the centre of band for low disorder (w) values, then it shifts 
towards the band edges for strong disorder. 


ie., number of nitrogen bases within a pitch of the helix, 
for various values of backbone disorder degree (w). It is 
clear that all the curves have the same general shape for 
the periodic as well as the random DNA sequences and 
the variation of 7 with v' is not monotonic. There exists 
a flat minima in these curves which indicates that at this 
point system is maximally extended. Now as we vary 
n (whatever be the disorder strength w is), 7 decreases. 


which indicates that system is less localized and effects of 
environmental fluctuations also becoming weaker. This 
behaviour can be explained easily, as we increase n we 
are allowing more channels for conduction between two 
adjacent pitches. As n increases, an electron can even¬ 
tually hops from one pitch to the next, galloping other 
nucleotides in that pitch. With increasing n, the length 
of this gallop also increases i.e., an electron gets the path 
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to bypass more number nucleotides as it move along the 
DNA chain. Because of this the effective length become 
shorter for an electron and it feels less disorder. Hence, 
first due to helical symmetry system become less local¬ 
ized and then due to conformation (n) it gets more and 
more extended. So, at this configuration system is hardly 
effected by external disturbances. This information can 
help to perform experiments on DNA in more easier way 
and reproducible results can be generated which is a chal¬ 
lenging task for a long time. 
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FIG. 5: (Color online). I-V response for two periodic se- 
qnences: poly(dA)-poly(dT) and poly(dG)-poly(dC) for five 
different disorder strength (w) at different values of n. For 
low disorder, cut off voltage reduces as we increase n, show¬ 
ing semiconducting behaviour. For strong disorder, current is 
considerably enhanced with increasing n giving a insulator to 
metallic transition. 

In Fig [5] we plot 7 vs. disorder strength (w), at a fixed 


value of t;'=0.3 eV, for several values of n. Here also as 
we increase n, 7 decreases for all values of w. But the 
variation of 7 with w is also not unidirectional. 7 reaches 
a peak value for disorder strength within 3<w<4 for all 
n values being considered. It signifies that at certain dis¬ 
order level, localization length becomes minimum, which 
implies that at this point system is most effected by exter¬ 
nal disturbances. This typical behaviour of localization 
is due to backbone structure of DNA . The effect of 
variation of n is also less for low disorder compared to 
higher ones. The effect of conformation (n) is maximum 
when the system is at its most localized state (3<w<4). 

In Fig [3] we show the variation of 7 with n. It also 
shows with increasing n, 7 decreases, except there is 
some different features around n=5 for one of the se¬ 
quence (ploy(dA)-poly(dT)). Though the reason is not 
yet clearly understood but it seems that system may has 
a critical configuration, at which it does feel environmen¬ 
tal effects most as we vary n. Because of that at n=5, 
7 increases instead of decreasing, showing it is the most 
localized configuration under appreciable disorder. This 
behaviour is not present in the other sequences, which 
shows different localization behaviour depending on se¬ 
quential variety. 

We also investigate localization behaviour with energy. 
In Fig|T]we plot the variation of 7 with energy for differ¬ 
ent values of n. The same thing is also happened here, 
with increasing n, 7 decreases. Though the rate of de¬ 
creasing is fast for small n (n=3, 5), then the variation 
of n become less effective for changes at higher values of 
n (n=7, 10). Effect is more prominent near centre of the 
bands for low disorder. As the disorder increases, effect 
of variation of n gradually delocalize towards the edges. 
At high disorder, variation is more sensitive around the 
edges of the band rather being at the centre. 

In Fig. [S] we plot I-V characteristics for the two peri¬ 
odic sequences for several values of n. We set the temper¬ 
ature at 0 K. To minimize the contact effects we choose 
tunnelling parameter r to be optimum i.e., x t 

between ds-DNA and the electrodes, where t is the hop¬ 
ping parameter for the electrodes [^. It is clear that 
effect of n is less at low disorder which is obvious be¬ 
cause at low disorder any path of charge conduction is 
equivalent as an electron feels almost no potential varia¬ 
tion. As the disorder increases effect of n becomes more 
distinctive. For strong disorder there is substantial vari¬ 
ation of potential at different sites and change in n gives 
an electron more number of shortcut pathways to move 
along the DNA chain. So, with increasing n, current is 
enhanced and the effect is sharp for high disorder values. 
For low disorder values cut-off voltage being reduced with 
increasing n, showing semiconductor-like transport. At 
high disorder for both the periodic sequences, current is 
considerably enhanced and almost linear response is ob¬ 
served at higher values of n, which indicates a transition 
from insulating to metallic phase. Our results are con¬ 
sistent with several experimental findings SilS- 
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IV. CONCLUDING REMARKS 

Till now different models have been used to study 
transport properties of DNA but none of these has taken 
into account of helical symmetry which is a basic feature 
of DNA structure. Using twisted ladder model we first in¬ 
corporate the helicity and then by varying the number of 
nucleotides within a pitch we try to model the conforma¬ 
tional variation of DNA. Though some calculations are 
present in the literature [26l - l^ but investigation within 
tight-binding framework is lacking. We report that de¬ 
pending on helical symmetry and conformation, localiza¬ 
tion properties can change considerably. The effect of 
conformation is less when environmental complications 
are small and increases with it. We have two interest¬ 
ing results. First one is by incorporating helical sym¬ 
metry and conformation we have been able to minimize 
the environmental effects to a great extent. It is clear 
from localization data that interplay of helical symme¬ 
try and conformation can provide some configurations 
where system is hardly disturbed by external agencies. 


If this information can be used correctly in experiments, 
we think the operation of such experiments would be¬ 
come less complicated. We investigated these properties 
in every aspect possible and it shows unambiguous vari¬ 
ation with conformational changes. The second result is, 
in presence of helical symmetry, depending on the coop¬ 
erative effect of backbone disorder and conformation sys¬ 
tem can undergo a transition from insulating to metallic 
phase as it is eminent from the TV responses of periodic 
sequences for higher disorder values. Whereas for low 
disorder with increasing n, cut-off voltage being reduced 
for semiconducting response. In summary, we can say 
that conformal changes have prominent effects on charge 
transport properties of DNA as it shows that DNA can 
be found in three different phases e.g., insulating, semi¬ 
conducting and metallic depending on the mutual varia¬ 
tion of environmental fluctuations and conformation. We 
hope in near future our results will be tested experimen¬ 
tally to find exact effects of helical symmetry as well as 
conformation on transport properties of DNA. 
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